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Abstract 

A search for the Standard Model Higgs boson in proton-proton collisions with the ATLAS detector 
at the LHC is presented. The datasets used correspond to integrated luminosities of approximately 
4.8 fb"' collected at ^ 7 TeV in 201 1 and 5.8 fb"^ at = 8 TeV in 201 2. Individual searches in the 
channels H-> ZZ^*^-> 4{, H-> yy and WW^*''^ evjiv in the 8 TeV data are combined with previously 
published results of searches for H^zz^*\ WW^*\ bb and t^t~ in the 7 TeV data and results from 
improved analyses of the zz^*'>^ Ai and yy channels in the 7 TeV data. Clear evidence for the 
production of a neutral boson with a measured mass of 126.0 + 0.4 (stat) + 0.4 (sys) GeV is presented. 
This observation, which has a significance of 5.9 standard deviations, corresponding to a background 
fluctuation probability of 1.7 x 10"^, is compatible with the production and decay of the Standard Model 
Higgs boson. 
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Abstract 

A search for the Standard Model Higgs boson in proton-proton collisions with the ATLAS detector at the 
LHC is presented. The datasets used correspond to integrated luminosities of approximately 4.8 fb ' collected at 
^/s = 7TeV in 2011 and 5.8 fb ' at = 8TeV in 2012. Individual searches in the channels H^ZZ^*^^4e, 
H^yy and //— > WW^*'— >evyuv in the 8TeV data are combined with previously published results of searches for 
H^ZZ^*\ hb and r+r" in the 7TeV data and results from improved analyses of the H^ZZ^*^^\l and 

yy channels in the 7 TeV data. Clear evidence for the production of a neutral boson with a measured mass of 
126.0 + 0.4 (stat) ± 0.4 (sys) GeV is presented. This observation, which has a significance of 5.9 standard devia- 
tions, corresponding to a background fluctuation probability of 1.7 x 10"'^, is compatible with the production and 
decay of the Standard Model Higgs boson. 



1. Introduction 

The Standard Model (SM) of particle physics |JJ 
has been tested by many experiments over the last four 
decades and has been shown to successfully describe 
high energy particle interactions. However, the mecha- 
nism that breaks electroweak symmetry in the SM has 
not been verified experimentally. This mechanism jsl- 
[loll . which gives mass to massive elementary particles, 
implies the existence of a scalar particle, the SM Higgs 
boson. The search for the Higgs boson, the only ele- 
mentary particle in the SM that has not yet been ob- 
served, is one of the highlights of the Large Hadron Col- 
lider 01 ill (LHC) physics programme. 

Indirect limits on the SM Higgs boson mass of rtiH < 
158 GeV at 95% confidence level (CL) have been set 
using global fits to precision electroweak results lll2ll . 
Direct searches at LEP IH], the Tevatron ifTj - fiell and 
the LHC HIHl have previously excluded, at 95% CL, 
a SM Higgs boson with mass below 600 GeV, apart from 
some mass regions between 1 16 GeV and 127 GeV. 

Both the ATLAS and CMS Collaborations reported 
excesses of events in their 2011 datasets of proton- 
proton (pp) collisions at centre-of-mass energy ^/s - 
7 TeV at the LHC, which were compatible with SM 
Higgs boson production and decay in the mass region 
124-126 GeV, with significances of2.9 and 3.1 standard 



deviations (cr), respectively The CDF and D0 

experiments at the Tevatron have also recently reported 
a broad excess in the mass region 120-1 35 GeV; using 
the existing LHC constraints, the observed local signifi- 
cances for oth - 125 GeV are 2.7 cr for CDF 1 141 , 1.1 o" 
for D0 iQ and 2.8 cr for flieir combination UM . 

The previous ATLAS searches in 4.6-4.8 fb ' of data 
at ^/s -1 TeV are combined here with new searches for 
ZZ^*'>^ 4/B yy and //-> WW**'^ evfiv in the 
5.8-5.9 fb"' of pp collision data taken at = 8 TeV 
between April and June 2012. 

The data were recorded with instantaneous luminosi- 
ties up to 6.8 X lO-'-' cm^^s"'; they are therefore af- 
fected by multiple pp collisions occurring in the same 
or neighbouring bunch crossings (pile-up). In the 7 TeV 
data, the average number of interactions per bunch 
crossing was approximately 10; the average increased 
to approximately 20 in the 8 TeV data. The reconstruc- 
tion, identification and isolation criteria used for elec- 
trons and photons in the 8 TeV data are improved, mak- 
ing the //— » ZZ'*^— » 4( and yy searches more ro- 
bust against the increased pile-up. These analyses were 
re-optimised with simulation and frozen before looking 
at the 8 TeV data. 



The symbol ( stands for electron or muon. 
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In the WW^*''^ (vtv channel, the increased pile- 
up deteriorates the event missing transverse momentum, 
£'™'*^ resolution, which results in significantly larger 
Drell-Yan background in the same-flavour final states. 
Since the eji channel provides most of the sensitivity 
of the search, only this final state is used in the anal- 
ysis of the 8 TeV data. The kinematic region in which 
a SM Higgs boson with a mass between 110 GeV and 
140 GeV is searched for was kept blinded during the 
analysis optimisation, until satisfactory agreement was 
found between the observed and predicted numbers of 
events in control samples dominated by the principal 
backgrounds. 

This Letter is organised as follows. The ATLAS de- 
tector is briefly described in Section |2l The simula- 
tion samples and the signal predictions are presented in 
Sectional The analyses of the ZZ**'^ 4^", yy 
and WW***^ ev^v channels are described in Sec- 
tions respectively. The statistical procedure used 
to analyse the results is summarised in Section |7] The 
systematic uncertainties which are correlated between 
datasets and search channels are described in Section[8j 
The results of the combination of all channels are re- 
ported in Section |9] while Section[TO]provides the con- 
clusions. 



2. The ATLAS detector 



The ATLAS detector \ is a multipurpose parti- 

cle physics apparatus with forward-backward symmet- 
ric cylindrical geometry. The inner tracking detector 
(ID) consists of a silicon pixel detector, a silicon mi- 
crostrip detector (SCT), and a straw-tube transition ra- 
diation tracker (TRT). The ID is surrounded by a thin su- 
perconducting solenoid which provides a 2 T magnetic 
field, and by high-granularity liquid-argon (LAr) sam- 
pling electromagnetic calorimetry. The electromagnetic 
calorimeter is divided into a central barrel (pseudora- 
pidit}0 It/I < 1.475) and end-cap regions on either end 
of the detector (1.375 < I/7I < 2.5 for the outer wheel 
and 2.5 < \rj\ < 3.2 for the inner wheel). In the region 
matched to the ID (I77I < 2.5), it is radially segmented 
into three layers. The first layer has a fine segmentation 



ATLAS uses a right-handed coordinate system with its origin at 
the nominal interaction point (IP) in the centre of the detector, and 
the z-axis along the beam line. The x-axis points from the IP to the 
centre of the LHC ring, and the y-axis points upwards. Cylindrical co- 
ordinates (r, 0) are used in the transverse plane, being the azimuthal 
angle around the beam line. Observables labelled "transverse" are 
projected into the x-y plane. The pseudorapidity is defined in terms 
of the polar angle 6 as = - In tan(6/2). 



in to facilitate ejy separation from and to improve 
the resolution of the shower position and direction mea- 
surements. In the region I77I < 1.8, the electromagnetic 
calorimeter is preceded by a presampler detector to cor- 
rect for upstream energy losses. An iron-scintillator/tile 
calorimeter gives hadronic coverage in the central ra- 
pidity range (I77I < 1.7), while a LAr hadronic end-cap 
calorimeter provides coverage over 1.5 < I77I < 3.2. The 
forward regions (3.2 < |?7| < 4.9) are instrumented with 
LAr calorimeters for both electromagnetic and hadronic 
measurements. The muon spectrometer (MS) surrounds 
the calorimeters and consists of three large air-core su- 
perconducting magnets providing a toroidal field, each 
with eight coils, a system of precision tracking cham- 
bers, and fast detectors for triggering. The combi- 
nation of all these systems provides charged particle 
measurements together with efficient and precise lepton 
and photon measurements in the pseudorapidity range 
I77I < 2.5. Jets and /i™"*^ are reconstructed using en- 
ergy deposits over the full coverage of the calorimeters, 
I77I < 4.9. 



3. Signal and background simulation samples 

The SM Higgs boson production processes con- 
sidered in this analysis are the dominant gluon fu- 
sion igg — > H, denoted ggF), vector-boson fusion 
iqq' qq'H, denoted VBF) and Higgs-strahlung 
(qq' WH,ZH, denoted WH/ZH). The small con- 
tribution from the associated production with a tl pair 
iqq/gg — > tlH, denoted ttH) is taken into account only 
in the //— > yy analysis. 

For the ggF process, the signal cross section is com- 
puted at up to next-to-next-to-leading order (NNLO) 
in QCD ll22U28l]. Next-to-leading order fNLO) elec- 
troweak (EW) corrections are applied ||29[ Isoll . as well 
as QCD soft-gluon re-summations at up to next-to-next- 



to-leading logarithm (NNLL) |31|. These calculations, 
which are described in Refs. ll32l - l35n . assume factori- 
sation between QCD and EW corrections. The trans- 
verse momentum, pT, spectrum of the Hig gs b oson in 
the ggF process follows the HqT calculation [36], which 
includes QCD corrections at NLO and QCD soft-gluon 
re-summations up to NNLL; the effects of finite quark 
masses are also taken into account 137 1. 

For the VBF process, full QCD and EW correc- 
tions up to NLO liiiSl] and approximate NNLO QCD 



corrections are used to calculate the cross sec- 

tion. Cross sections of the associated WH/ZH processes 
(VH) are calculated including QCD corrections up to 
NNLO liliQ and EW corrections up to NLO S. 
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The cross sections for the ttH process are estimated up 
toNLO QCD 147^5111. 

The total cross sections for SM Higgs boson produc- 
tion at the LHC with nin - 125 GeV are predicted to 
be 17.5 pb for ^fs - 7TeV and 22.3 pb for ^fs = 
STeViillll. 

The branching ratios of the SM Higgs boson as a 
function of mh, as well as their uncertainties, are calcu- 
lated using the HDECAY [54] and PROPHECY4F [55, 
l56ll programs and are taken from Refs. |52, 53 1. The 
interference in the ZZ*^*^^ At final states with iden- 
tical leptons is taken into account Ii53l,55.,56.] . 



Table 1 : Event generators used to model the signal and background 
processes. "PYTHIA" indicates that PYTHIA6 and PYTHIAS are 
used for simulations of ■\/s = 1 TeV and = 8 TeV data, respec- 
tively. 



Process 


Generator 


ggF, VBF 


POWHEG [57, 58J+PYTHIA 


WH, ZH, ttH 


PYTHIA 


W+jets, Z/y* +jets 


ALPGEN [591+HERWIG 


t't, tW, tb 


MC@NLO [6gi+HERWIG 


tqb 


AcerMC [6r|+PYTHIA 


qq WW 


MC@NLO+HERWIG 


gg^WW 


gg2WW [62]+HERWIG 


qq ZZ 


POWHEG [631+PYTHIA 


gg^zz 


gg2ZZ [64]+HERWIG 


wz 


MadGraph+PYTHIA, HERWIG 


Wy+jets 


ALPGEN+HERWIG 


Wy' [65] 


MadGraph+PYTHIA 


qqigg rr 


SHERPA 



The event generators used to model signal and back- 
ground processes in samples of Monte Carlo (MC) sim- 
ulated events are listed in Table [T] The normalisations 
of the generated samples are obtained from the state of 
the art calculations described above. Several different 
programs are used to generate the hard-scattering pro- 
cesses. To generate parton showers and their hadroni- 
sation, and to^simulate the underlying event 
PYTHIA6 169] (for 7 TeV samples and 8 TeV sam- 
ples produced with MadGraph \JQ^ 21} or AcerMC) or 
PYTHIAS 111 (for other 8 TeV samples) are used. Al- 
ternatively, HERWIG fr^ or SHERPA fT?] are used 
to generate and hadronise parton showers, with the 
HERWIG underlying event simulation performed using 
JIMMY f7?]. When PYTHIA6 or HERWIG are used, 
TAUOLA [76J and PHOTOS |77] are employed to de- 
scribe tau lepton decays and additional photon radiation 
from charged leptons, respectively. 



The following parton distribution function (PDF) sets 
are used: CTIO |78] for the POWHEG, MC@NLO, 
gg2WW and gg2ZZ samples; CTEQ6L1 S for the 
PYTHIA8, ALPGEN, AcerMC, MadGraph, HERWIG 
and SHERPA samples; and MRSTMCal |80] for the 
PYTHIA6 samples. 

Acceptances and efficiencies are obtained mostly 
from full simulations of the ATLAS detector Ii8 LI us- 
ing Geant4 11821] . These simulations include a realistic 
modelling of the pile-up conditions observed in the data. 
Corrections obtained from measurements in data are ap- 
plied to account for small differences between data and 
simulation (e.g. large samples of W, Z and decays 
are used to derive scale factors for lepton reconstruction 
and identification efficiencies). 

A. H ^ ZZ^*^ ^ 4f channel 

The search for the SM Higgs boson through the 
decay H — > ZZ**' — > At, where { - e or ju, pro- 
vides good sensitivity over a wide mass range (110- 
600 GeV), largely due to the excellent momentum reso- 
lution of the ATLAS detector. This analysis searches 
for Higgs boson candidates by selecting two paks of 
isolated leptons, each of which is comprised of two lep- 
tons with the same flavour and opposite charge. The 
expected cross section times branching ratio for the pro- 
cess H ZZ^*> At with niH = 125 GeV is 2.2 fb for 
= 7 TeV and 2.8 fb for V* = 8 TeV. 

The largest background comes from continuum 
(Z^*^ ly*){Z^*^ ly*) production, referred to hereafter as 
ZZ'*'. For low masses there are also important back- 
ground contributions from Z + jets and tt production, 
where charged lepton candidates arise either from de- 
cays of hadrons with b- or c-quark content or from mis- 
identification of jets. 

The 7 TeV data have been re-analysed and combined 
with the 8 TeV data. The analysis is improved in several 
aspects with respect to Ref. [83,1 to enhance the sensitiv- 
ity to a low-mass Higgs boson. In particular, the kine- 
matic selections are revised, and the 8 TeV data anal- 
ysis benefits from improvements in the electron recon- 
struction and identification. The expected signal sig- 
nificances for a Higgs boson with mn - 125 GeV are 
1.6 cr for the 7 TeV data (to be compared with 1.25 cr 
in Ref. HI) and 2. 1 cr for the 8 TeV data. 



4.1. Event selection 

The data are selected using single-lepton or dilepton 
triggers. For the single-muon trigger, the pi threshold 
is 18 GeV for the 7 TeV data and 24 GeV for the 8 TeV 
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data, while for the single-electron trigger the transverse 
energy, Ej, threshold varies from 20 GeV to 22 GeV for 
the 7 TeV data and is 24 GeV for the 8 TeV data. For 
the dielectron triggers, the thresholds are 12 GeV for 
both electrons. For the dimuon triggers, the thresholds 
for the 7 TeV data are 10 GeV for each muon, while 
for the 8 TeV data the thresholds are 13 GeV. An addi- 
tional asymmetric dimuon trigger is used in the 8 TeV 
data with thresholds 18 GeV and 8 GeV for the leading 
and sub-leading muon, respectively. 

Muon candidates are formed by matching recon- 
structed ID tracks with either a complete track or a 
track-segment reconstructed in the MS ll84ll . The muon 
acceptance is extended with respect to Ref. fs?] us- 
ing tracks reconstructed in the forward region of the 
MS (2.5 < It/I < 2.7), which is outside the ID cov- 
erage. If both an ID and a complete MS track are 
present, the two independent momentum measurements 
are combined; otherwise the information of the ID or 
the MS is used alone. Electron candidates must have a 
well-reconstructed ID track pointing to an electromag- 
netic calorimeter cluster and the cluster should satisfy 
a set of identification criteria fsS*] that require the lon- 
gitudinal and transverse shower profiles to be consis- 
tent with those expected for electromagnetic showers. 
Tracks associated with electromagnetic clusters are fit- 
ted using a Gaussian-Sum Filter |86], which allows for 
bremsstrahlung energy losses to be taken into account. 

Each electron (muon) must satisfy pj > 7 GeV 
(Pt > 6 GeV) and be measured in the pseudorapidity 
range It/I < 2.47 (|t/| < 2.7). All possible quadruplet 
combinations with same-flavour opposite-charge lepton 
pairs are then formed. The most energetic lepton in the 
quadruplet must satisfy pj > 20 GeV, and the second 
(third) lepton in pj order must satisfy pj > 15 GeV 
(Pt > 10 GeV). At least one of the leptons must 
satisfy the single-lepton trigger or one pair must sat- 
isfy the dilepton trigger requirements. The leptons are 
required to be separated from each other by AR - 
yjiAriY + (A0)2 > 0.1 if they are of the same flavour 
and by A/? > 0.2 otherwise. The longitudinal impact pa- 
rameters of the leptons along the beam axis are required 
to be within 10 mm of the reconstructed primary vertex. 
The primary vertex used for the event is defined as the 
reconstructed vertex with the highest 2 p\ of associated 
tracks and is required to have at least three tracks with 
Pi > 0.4 GeV. To reject cosmic rays, muon tracks are 
required to have a transverse impact parameter, defined 
as the distance of closest approach to the primary vertex 
in the transverse plane, of less than 1 mm. 

The same-flavour and opposite-charge lepton pair 



with an invariant mass closest to the Z boson mass {niz) 
in the quadruplet is referred to as the leading lepton pair 
Its invariant mass, denoted by nin, is required to be 
between 50 GeV and 106 GeV. The remaining same- 
flavour, opposite-charge lepton pair is the sub-leading 
lepton pair. Its invariant mass, 0134, is required to be 
in the range Wmin < '«34 < 115 GeV, where the value 
of OTmin depends on the reconstructed four-lepton in- 
variant mass, mm. The value of mmi„ varies monoton- 
ically from 17.5 GeV at = 120 GeV to 50 GeV 
at - 190 GeV IstIi and is constant above this 
value. All possible lepton pairs in the quadruplet that 
have the same flavour and opposite charge must satisfy 
mtt > 5 GeV in order to reject backgrounds involving 
the production and decay of mesons. If two or more 
quadruplets satisfy the above selection, the one with the 
highest value of OT34 is selected. Four different analysis 
sub-channels, 4e, 2e2yu, 2jJ2e and 4jU, arranged by the 
flavour of the leading lepton pair, are defined. 

Non-prompt leptons from heavy flavour decays, elec- 
trons from photon conversions and jets mis-identified 
as electrons have broader transverse impact parameter 
distributions than prompt leptons from Z boson decays 
and/or are non-isolated. Thus, the Z+jets and tt back- 
ground contributions are reduced by applying a cut on 
the transverse impact parameter significance, defined as 
the transverse impact parameter divided by its uncer- 
tainty, do/cTd^. This is required to be less than 3.5 (6.5) 
for muons (electrons). The electron impact parameter is 
aff'ected by bremsstrahlung and thus has a broader dis- 
tribution. 

In addition, leptons must satisfy isolation require- 
ments based on tracking and calorimetric information. 
The normalised track isolation discriminant is defined 
as the sum of the transverse momenta of tracks inside a 
cone of size AR = 0.2 around the lepton direction, ex- 
cluding the lepton track, divided by the lepton pj. The 
tracks considered in the sum are those compatible with 
the lepton vertex and have px > 0.4 GeV (pj > 1 GeV) 
in the case of electron (muon) candidates. Each lepton 
is required to have a normalised track isolation smaller 
than 0.15. The normalised calorimetric isolation for 
electrons is computed as the sum of the Ej of positive- 
energy topological clusters L88J with a reconstructed 
barycentre falling within a cone of size AR = 0.2 around 
the candidate electron cluster, divided by the electron 
Ej. The algorithm for topological clustering suppresses 
noise by keeping ceUs with a significant energy deposit 
and their neighbours. The summed energy of the cells 
assigned to the electron cluster is excluded, while a cor- 
rection is applied to account for the electron energy de- 
posited outside the cluster The ambient energy deposi- 
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tion in the event from pile-up and the underlying event 
is accounted for using a calculation of the median trans- 
verse energy density from 1ow-/9t jets li^S^l- The nor- 
malised calorimetric isolation for electrons is required 
to be less than 0.20. The normalised calorimetric isola- 
tion discriminant for muons is defined by the ratio to the 
Pt: of the muon of the Ej sum of the calorimeter cells 
inside a cone of size AR = 0.2 around the muon direc- 
tion minus the energy deposited by the muon. Muons 
are required to have a normalised calorimetric isolation 
less than 0.30 (0. 15 for muons without an associated ID 
track). For both the track- and calorimeter-based isola- 
tion, any contributions arising from other leptons of the 
quadruplet are subtracted. 

The combined signal reconstruction and selection ef- 
ficiencies for a SM Higgs with niH = 125 GeV for the 
7 TeV (8 TeV) data are 37% (36%) for the V channel, 
20% (22%) for the leliijlfjle channels and 15% (20%) 
for the 4e channel. 

The A{ invariant mass resolution is improved by ap- 
plying a Z-mass constrained kinematic fit to the leading 
lepton pair for m^f < 190 GeV and to both lepton pairs 
for higher masses. The expected width of the recon- 
structed mass distribution is dominated by the experi- 
mental resolution for mn < 350 GeV, and by the natu- 
ral width of the Higgs boson for higher masses (30 GeV 
at niH - 400 GeV). The typical mass resolutions for 
niH = 125 GeV ai-e 1.7 GeV, 1.7 GeV/2.2 GeV and 
2.3 GeV for the 4yU, 2e2/i/2yu2e and 4e sub-channels, re- 
spectively. 

4.2. Background estimation 

The expected background yield and composition are 
estimated using the MC simulation normalised to the 
theoretical cross section for ZZ*^*^ production and by 
methods using control regions from data for the Z -H jets 
and ff processes. Since the background composition de- 
pends on the flavour of the sub-leading lepton pair, dif- 
ferent approaches are taken for the ii+^fj. and the tt+ee 
final states. The transfer factors needed to extrapolate 
the background yields from the control regions defined 
below to the signal region are obtained from the MC 
simulation. The MC description of the selection effi- 
ciencies for the diff'erent background components has 
been verified with data. 

The reducible II + jjjj background is dominated by ff 
and Z -H jets (mostly Zbb) events. A control region is de- 
fined by removing the isolation requirement on the lep- 
tons in the sub-leading pair, and by requiring that at least 
one of the sub-leading muons fails the transverse impact 
parameter significance selection. These modifications 
remove ZZ*^** contributions, and allow both the ff and 



Z H- jets backgrounds to be estimated simultaneously us- 
ing a fit to the mi2 distribution. The ff background con- 
tribution is cross-checked by selecting a control sample 
of events with an opposite charge efi pair with an invari- 
ant mass between 50 GeV and 106 GeV, accompanied 
by an opposite-charge muon pair Events with a Z can- 
didate decaying to a pair of electrons or muons in the 
aforementioned mass range are excluded. Isolation and 
transverse impact parameter significance requirements 
are applied only to the leptons of the eju pair. 

In order to estimate the reducible {{ + ee background, 
a control region is formed by relaxing the selection cri- 
teria for the electrons of the sub-leading pair. The diff'er- 
ent sources of electron background are then separated 
into categories consisting of non-prompt leptons from 
heavy flavour decays, electrons from photon conver- 
sions and jets mis-identified as electrons, using appro- 
priate discriminating variables 119 ill . This method allows 
the sum of the Z H- jets and ff background contributions 
to be estimated. As a cross-check, the same method is 
also applied to a similar control region containing same- 
charge sub-leading electron pairs. An additional cross- 
check of the ll+ee background estimation is performed 
by using a control region with same-charge sub-leading 
electron pairs, where the three highest pi leptons satisfy 
all the analysis criteria whereas the selection cuts are re- 
laxed for the remaining electrons. All the cross-checks 
yield consistent results. 



Table 2: Summary of the estimated numbers of Z + jets and tt back- 
ground events, for the Vs = 7 TeV and = 8 TeV data in the entire 
phase-space of the analysis after the kinematic selections described in 
the text. The backgrounds are combined for the 2//2e and 4e channels, 
as discussed in the text. The first uncertainty is statistical, while the 
second is systematic. 



Background 


Estimated 




numbers of events 




VI = 7 TeV = 8 TeV 




Z+jets 


0.3+0.1 ±0.1 0.5+0.1+0.2 


tt 


0.02+0.02+0.01 0.04+0.02+0.02 


2e2yU 


Z+jets 


0.2+0.1+0.1 0.4+0.1+0.1 


ti 


0.02+0.01+0.01 0.04+0.01+0.01 


lya.e 


Z+jets, tt 


2.6+0.4+04 4.9+0.8+0.7 


4e 


Z+jets, tt 


3.1+0.6+0.5 3.9+0.7+0.8 



The data-driven background estimates are sum- 
marised in Table |2] The distribution of m34, for events 
selected by the analysis except that the isolation and 
transverse impact parameter requirements for the sub- 
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leading lepton pair are removed, is presented in Fig.[T] 
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Figure 1: Invariant mass distribution of the sub-leading lepton pair 
(m34) for a sample defined by the presence of a Z boson candidate and 
an additional same-flavour electron or muon pair, for the coinbination 
of = 7 TeV and = 8 TeV data in the entire phase-space of the 
analysis after the kinematic selections described in the text. Isolation 
and transverse impact parameter significance requirements are applied 
to the leading lepton pair only. The MC is normalised to the data- 
diiven background estiinations. The relatively small contribution of a 
SM Higgs with m// = 125 GeV in this sample is also shown. 



4.3. Systematic uncertainties 

The uncertainties on the integrated luminosities are 
determined to be 1.8% for the 7 TeV data and 3.6% 
for the 8 TeV data using the techniques described in 
Ref. {9^]. 

The uncertainties on the lepton reconstruction and 
identification efficiencies and on the momentum scale 
and resolution are determined using samples of W, 
Z and y/i/r decays 1184 18511 . The relative uncertainty 
on the signal acceptance due to the uncertainty on 
the muon reconstruction and identification efficiency is 
+0.7% (+0.5%/±0.5%) for the Aji {leljillixle) chan- 
nel for m/ii - 600 GeV and increases to +0.9% 
(+0.8%/±0.5%) for = 115 GeV. Similarly, the 
relative uncertainty on the signal acceptance due to the 
uncertainty on the electron reconstruction and identifi- 
cation efficiency is +2.6% (+1.7%/+ 1.8%) for the 4e 
(2e2///2//2e) channel for mn( - 600 GeV and reaches 
±8.0% (+2.3%/±7.6%) for niv = 115 GeV. The un- 
certainty on the electron energy scale results in an un- 
certainty of +0.7% (+0.5%/+0.2%) on the mass scale 
of the m-At distribution for the 4e (2e2fi/2fi2e) channel. 
The impact of the uncertainties on the electron energy 



resolution and on the muon momentum resolution and 
scale are found to be negligible. 

The theoretical uncertainties associated with the sig- 
nal are described in detail in Section [8] For the SM 
ZZ*^*' background, which is estimated from MC simula- 
tion, the uncertainty on the total yield due to the QCD 
scale uncertainty is ±5%, while the effect of the PDF 
and as uncertainties is ±4% (+8%) for processes initi- 
ated by quarks (gluons) 115 311 . In addition, the depen- 
dence of these uncertainties on the four-lepton invariant 
mass spectrum has been taken into account as discussed 
in Ref. 115 311 . Though a small excess of events is ob- 
served for m4i > 160 GeV, the measured ZZ*^*' — > 4i 
cross section [93] is consistent with the SM theoreti- 
cal prediction. The impact of not using the theoretical 
constraints on the ZZ**^ yield on the search for a Higgs 
boson with niH < 2mz has been studied in Ref. Ii87il and 
has been found to be negligible . The impact of the in- 
terference between a Higgs signal and the non-resonant 
gg — > ZZ^*' background is small and becomes negligi- 
ble for niH < 2mz 1.94.1 . 
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Figure 2: The distribution of the four-lepton invariant mass, m^i, for 
the selected candidates, compared to the background expectation in 
the 80-250 GeV mass range, for the combination of the = 7 TeV 
and = 8 TeV data. The signal expectation for a SM Higgs with 
mn = 125 GeV is also shown. 



4.4. Results 

The expected distributions of m4( for the background 
and for a Higgs boson signal with nin = 125 GeV are 
compared to the data in Fig. |2] The numbers of ob- 
served and expected events in a window of ±5 GeV 
around mn - 125 GeV are presented for the combined 
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7 TeV and 8 TeV data in Table [3] The distribution of 
the versus mi2 invariant mass is shown in Fig. [3] 
The statistical interpretation of the excess of events near 
m4c = 125 GeV in Fig.|2]is presented in Section|9] 



Table 3: The numbers of expected signal (niM = 125 GeV) and back- 
ground events, together with the numbers of observed events in the 
data, in a window of size ±5 GeV around 125 GeV, for the combined 
V? = 7 TeV and V* = 8 TeV data. 





Signal 


ZZM 


Z + jets, tt 


Observed 




2.09±0.30 


1.12±0.05 


0.13±0.04 


6 


lelfi/lfile 


2.29± 0.33 


0.80±0.05 


1.27±0.19 


5 


4e 
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0.44±0.04 


1.09±0.20 
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Figure 3: Distribution of the 0134 versus the mi2 invariant mass, be- 
fore the application of the Z-mass constrained kinematic fit, for the 
selected candidates in the m4( range 120-130 GeV. The expected 
distributions for a SM Higgs with mu = 125 GeV (the sizes of the 
boxes indicate the relative density) and for the total background (the 
intensity of the shading indicates the relative density) are also shown. 



5. //-» 77 channel 

The search for the SM Higgs boson through the de- 
cay //— >77 is performed in the mass range between 
110 GeV and 150 GeV. The dominant background is 
SM diphoton production (yy); contributions also come 
from 7+jet and jet-t-jet production with one or two jets 
mis-identified as photons (y j and jj) and from the Drell- 
Yan process. The 7 TeV data have been re-analysed and 
the results combined with those from the 8 TeV data. 
Among other changes to the analysis, a new category 
of events with two jets is introduced, which enhances 



the sensitivity to the VBF process. Higgs boson events 
produced by the VBF process have two forward jets, 
originating from the two scattered quarks, and tend to 
be devoid of jets in the central region. Overall, the sen- 
sitivity of the analysis has been improved by about 20% 
with respect to that described in Ref. 1.95,1 . 

5.7. Event selection 

The data used in this channel are selected using 
a diphoton trigger i96ll . which requires two clusters 
formed from energy depositions in the electromagnetic 
calorimeter An Ej threshold of 20 GeV is applied to 
each cluster for the 7 TeV data, while for the 8 TeV 
data the thresholds are increased to 35 GeV on the lead- 
ing (the highest Ej) cluster and to 25 GeV on the sub- 
leading (the next-highest Ej) cluster In addition, loose 
criteria are applied to the shapes of the clusters to match 
the expectations for electromagnetic showers initiated 
by photons. The efficiency of the trigger is greater than 
99% for events passing the final event selection. 

Events are required to contain at least one recon- 
structed vertex with at least two associated tracks with 
Pj > 0.4 GeV, as well as two photon candidates. Pho- 
ton candidates are reconstructed in the fiducial region 
\ri\ < 2.37, excluding the calorimeter barrel/end-cap 
transition region 1.37 < I?;! < 1.52. Photons that convert 
to electron-positron pairs in the ID material can have 
one or two reconstructed tracks matched to the clusters 
in the calorimeter The photon reconstruction efficiency 
is about 97% for Ej > 30 GeV. 

In order to account for energy losses upstream of the 
calorimeter and energy leakage outside of the cluster, 
MC simulation results are used to calibrate the energies 
of the photon candidates; there are separate calibrations 
for unconverted and converted candidates. The calibra- 
tion is refined by applying ?7-dependent correction fac- 
tors, which are of the order of ± 1 %, determined from 
measured Z— > e^e^ events. The leading (sub-leading) 
photon candidate is required to have Ej > 40 GeV 
(30 GeV). 

Photon candidates are required to pass identification 
criteria based on shower shapes in the electromagnetic 
calorimeter and on energy leakage into the hadronic 
calorimeter 19711 . For the 7 TeV data, this information is 
combined in a neural network, tuned to achieve a sim- 
ilar jet rejection as the cut-based selection described in 
Ref. 195], but with higher photon efficiency. For the 
8 TeV data, cut-based criteria are used to ensure reliable 
photon performance for recently-recorded data. This 
cut-based selection has been tuned to be robust against 
pile-up by relaxing requirements on shower shape cri- 
teria more susceptible to pile-up, and tightening others. 
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The photon identification efficiencies, averaged over 
range from 85% to above 95% for the Ej range under 
consideration. 

To further suppress the jet background, an isolation 
requirement is applied. The isolation transverse en- 
ergy is defined as the sum of the transverse energy 
of positive-energy topological clusters, as described in 
Section m within a cone of size AR = 0.4 around the 
photon candidate, excluding the region within 0.125 x 
0.175 in A77X around the photon barycentre. The dis- 
tributions of the isolation transverse energy in data and 
simulation have been found to be in good agreement us- 
ing electrons from Z— > e^e" events and photons from 
Z l^Vj events. Remaining small differences are 
taken into account as a systematic uncertainty. Photon 
candidates are required to have an isolation transverse 
energy of less than 4 GeV. 

5.2. Invariant mass reconstruction 

The invariant mass of the two photons is evaluated us- 
ing the photon energies measured in the calorimeter, the 
azimuthal angle ^ between the photons as determined 
from the positions of the photons in the calorimeter, and 
the values of r\ calculated from the position of the identi- 
fied primary vertex and the impact points of the photons 
in the calorimeter 

The primary vertex of the hard interaction is identi- 
fied by combining the following information in a global 
likelihood: the directions of flight of the photons as 
determined using the longitudinal segmentation of the 
electromagnetic calorimeter (calorimeter pointing), the 
parameters of the beam spot, and the 2 p\ of the tracks 
associated with each reconstructed vertex. In addition, 
for the 7 TeV data analysis, the reconstructed conver- 
sion vertex is used in the likelihood for converted pho- 
tons with tracks containing hits in the silicon layers of 
the ID. The calorimeter pointing is sufficient to ensure 
that the contribution of the opening angle between the 
photons to the mass resolution is negligible. Using the 
calorimeter pointing alone, the resolution of the vertex z 
coordinate is ~ 15 mm, improving to ~ 6 mm for events 
with two reconstructed converted photons. The tracking 
information from the ID improves the identification of 
the vertex of the hard interaction, which is needed for 
the jet selection in the 2-jet category. 

With the selection described in Section 15.11 in the 
diphoton invariant mass range between 100 GeV and 
160 GeV, 23788 and 35251 diphoton candidates are ob- 
served in the 7 TeV and 8 TeV data samples, respec- 
tively. 

Data-driven techniques llsisll are used to estimate the 
numbers of 77, yj and jj events in the selected sam- 



ple. The contribution from the Drell-Yan background 
is determined from a sample of Z— > e^e' decays in data 
where either one or both electrons pass the photon selec- 
tion. The measured composition of the selected sample 
is approximately 74%, 22%, 3% and 1% for the yy, yj, 
jj and Drell-Yan processes, respectively, demonstrating 
the dominance of the irreducible diphoton production. 
This decomposition is not directly used in the signal 
search; however, it is used to study the parameterisation 
of the background modelling. 

5.3. Event categorisation 

To increase the sensitivity to a Higgs boson signal, 
the events are separated into ten mutually exclusive cat- 
egories having different mass resolutions and signal-to- 
background ratios. An exclusive category of events con- 
taining two jets improves the sensitivity to VBF. The 
other nine categories are defined by the presence or not 
of converted photons, of the selected photons, and pjt, 
the componenji of the diphoton pj that is orthogonal to 
the axis defined b y the difference between the two pho- 
ton momenta Slloil - 



Jets are reconstructed OlOlll using the anti-^, algo- 
rithm 1 102] with radius parameter R - 0.4. At least 
two jets with \ri\ < 4.5 and pj > 25 GeV are re- 
quired in the 2-jet selection. In the analysis of the 8 TeV 
data, the pj threshold is raised to 30 GeV for jets with 
2.5 < I77I < 4.5. For jets in the ID acceptance (I77I < 2.5), 
the fraction of the sum of the pj of tracks, associated 
with the jet and matched to the selected primary vertex, 
with respect to the sum of the p-Y of tracks associated 
with the jet (jet vertex fraction, JVF) is required to be 
at least 0.75. This requirement on the JVF reduces the 
number of jets from proton-proton interactions not asso- 
ciated with the primary vertex. Motivated by the VBF 
topology, three additional cuts are applied in the 2-jet 
selection: the difference of the pseudorapidity between 
the leading and sub-leading jets (tag jets) is required to 
be larger than 2.8, the invariant mass of the tag jets has 
to be larger than 400 GeV, and the azimuthal angle dif- 
ference between the diphoton system and the system of 
the tag jets has to be larger than 2.6. About 70% of the 
signal events in the 2-jet category come from the VBF 
process. 

The other nine categories are defined as follows: 
events with two unconverted photons are separated into 
unconverted central (\ri\ < 0.75 for both candidates) and 
unconverted rest (all other events), events with at least 



are the transverse momenta of the two photons. 



p^- 1, where p^' and p^- 
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Table 4: Number of events in the data (Nd) and expected number of 
signal events (Ws) for nifj = 126.5 GeV from the //— > yy analysis, for 
each category in the mass range 100-160 GeV. The mass resolution 
FWHM (see text) is also given for the STeVdata. The Higgs boson 
production cross section multiplied by the branching ratio into two 
photons ((TXB(H — > yy)) is listed for inn = 126.5 GeV. The statistical 
uncertainties on Ns and FWHM are less than 1 %. 





7TeV 


8TeV 


a- X B(H -> yy) [lb] 


39 


50 


FWHM 


Category 


Nd 


Ns 


Nd 


Ns 


[GeV] 


Unconv. central, low pxt 


2054 


10.5 


2945 


14.2 


3.4 


Unconv. central, high pxt 


97 


1.5 


173 


2.5 


3.2 


Unconv. rest, low pxt 


7129 


21.6 


12136 


30.9 


3.7 


Unconv. rest, high pxt 


444 


2.8 


785 


5.2 


3.6 


Conv. central, low /?Tt 


1493 


6.7 


2015 


8.9 


3.9 


Conv. central, high pj^ 


77 


1.0 


113 


1.6 


3.5 


Conv. rest, low pji 


8313 


21.1 


11099 


26.9 


4.5 


Conv. rest, high /?Tt 


501 


2.7 


706 


4.5 


3.9 


Conv. transition 


3591 


9.5 


5140 


12.8 


6.1 


2-jet 


89 


2.2 


139 


3.0 


3.7 


All categories (inclusive) 


23788 


79.6 


35251 


110.5 


3.9 



one converted photon are separated into converted cen- 
tral {\r]\ < 0.75 for both candidates), converted transi- 
tion (at least one photon with 1.3 < I77I < 1.75) and 
converted rest (all other events). Except for the con- 
verted transition category, each category is further di- 
vided by a cut at pn- 60 GeV into two categories, low 
pn and high pn. MC studies show that signal events, 
particularly those produced via VBF or associated pro- 
duction {WHjZH and ttH), have on average larger pn 
than background events. The number of data events in 
each category, as well as the sum of all the categories, 
which is denoted inclusive, are given in Table H) 

5.4. Signal modelling 

The description of the Higgs boson signal is obtained 
from MC, as described in Section [3] The cross sections 
multiplied by the branching ratio into two photons are 
given in Table |4] for nin - 126.5 GeV. The number of 
signal events produced via the ggF process is rescaled 
to take into account the expected destructive interfer- 
ence between the gg — > yy continuum background and 
ggF II103II . leading to a reduction of the production rate 
by 2-5% depending on nin and the event category. For 
both the 7 TeV and 8 TeV MC samples, the fractions of 
ggF, VBF, WH, ZH and ttH production are approxi- 
mately 88%, 7%, 3%, 2% and 0.5%, respectively, for 
niH = 126.5 GeV. 

In the simulation, the shower shape distributions 
are shifted slightly to improve the agreement with the 
data II97II . and the photon energy resolution is broad- 
ened (by approximately 1% in the barrel calorimeter 



and 1.2-2.1% in the end-cap regions) to account for 
small differences observed between e^e" data and 
MC events. The signal yields expected for the 7 TeV 
and 8 TeV data samples are given in Table H) The over- 
all selection efficiency is about 40%. 

The shape of the invariant mass of the signal in each 
category is modelled by the sum of a Crystal Ball func- 
tion 1.104.1 . describing the core of the distribution with 
a width (TcB, and a Gaussian contribution describing 
the tails (amounting to <10%) of the mass distribution. 
The expected full-width-at-half-maximum (FWHM) is 
3.9 GeV and ctcb is 1.6 GeV for the inclusive sample. 
The resolution varies with event category (see Table|4|; 
the FWHM is typically a factor 2.3 larger than ctcb- 

5.5. Background modelling 

The background in each category is estimated from 
data by fitting the diphoton mass spectrum in the mass 
range 100-160 GeV with a selected model with free pa- 
rameters of shape and normalisation. Different models 
are chosen for the different categories to achieve a good 
compromise between limiting the size of a potential bias 
while retaining good statistical po\y er. A fourth-order 
Bernstein polynomial function 1 105] is used for the un- 
converted rest (low Pti), converted rest (low pn) and in- 
clusive categories, an exponential function of a second- 
order polynomial for the unconverted central (low pjt), 
converted central (low pxt) and converted transition cat- 
egories, and an exponential function for all others. 

Studies to determine the potential bias have been per- 
formed using large samples of simulated background 
events complemented by data-driven estimates. The 
background shapes in the simulation have been cross- 
checked using data from control regions. The poten- 
tial bias for a given model is estimated, separately for 
each category, by performing a maximum likelihood fit 
to large samples of simulated background events in the 
mass range 100-160 GeV, of the sum of a signal plus 
the given background model. The signal shape is taken 
to follow the expectation for a SM Higgs boson; the sig- 
nal yield is a free parameter of the fit. The potential bias 
is defined by the largest absolute signal yield obtained 
from the likelihood fit to the simulated background sam- 
ples for hypothesised Higgs boson masses in the range 
1 10-150 GeV. A pre-selection of background parame- 
terisations is made by requiring that the potential bias, 
as defined above, is less than 20% of the statistical un- 
certainty on the fitted signal yield. The pre-selected pa- 
rameterisation in each category with the best expected 
sensitivity for niH = 125 GeV is selected as the back- 
ground model. 
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The largest absolute signal yield as defined above is 
taken as the systematic uncertainty on the background 
model. It amounts to ±(0.2-4.6) and +(0.3-6.8) events, 
depending on the category for the 7 TeV and 8 TeV data 
samples, respectively. In the final fit to the data (see 
Section 15^ a signal-like term is included in the likeli- 
hood function for each category. This term incorporates 
the estimated potential bias, thus providing a conserva- 
tive estimate of the uncertainty due to the background 
modelling. 

5.6. Systematic uncertainties 

Hereafter, in cases where two uncertainties are 
quoted, they refer to the 7 TeV and 8 TeV data, respec- 
tively. The dominant experimental uncertainty on the 
signal yield (±8%, +11%) comes from the photon re- 
construction and identification efliciency, which is es- 
timated with data using electrons from Z decays and 
photons from Z £^i^y events. Pile-up modelling 
also affects the expected yields and contributes to the 
uncertainty (+4%). Further uncertainties on the sig- 
nal yield are related to the trigger (± 1 %), photon isola- 
tion (±0.4%, ±0.5%) and luminosity (±1.8%, ±3.6%). 
Uncertainties due to the modelling of the underlying 
event are ±6% for VBF and ±30% for other produc- 
tion processes in the 2-jet category. Uncertainties on the 
predicted cross sections and branching ratio are sum- 
marised in Section|8] 

The uncertainty on the expected fractions of signal 
events in each category is described in the following. 
The uncertainty on the knowledge of the material in 
front of the calorimeter is used to derive the amount of 
possible event migration between the converted and un- 
converted categories (±4%). The uncertainty from pile- 
up on the population of the converted and unconverted 
categories is ±2%. The uncertainty from the jet energy 
scale (JES) amounts to up to ±19% for the 2-jet cate- 
gory, and up to ±4% for the other categories. Uncertain- 
ties from the JVF modelling are ±12% (for the 8 TeV 
data) for the 2-jet category, estimated from Z+2-jets 
events by comparing data and MC. Different PDFs and 
scale variations in the HqT calculations are used to de- 
rive possible event migration among categories (±9%) 
due to the modelling of the Higgs boson kinematics. 

The total uncertainty on the mass resolution is ±14%. 
The dominant contribution (±12%) comes from the un- 
certainty on the energy resolution of the calorimeter, 
which is determined from Z— > e^e^ events. Smaller 
contributions come from the imperfect knowledge of the 
material in front of the calorimeter, which affects the ex- 
trapolation of the calibration from electrons to photons 
(±6%), and from pile-up (±4%). 
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Figure 4: The distributions of tlie invariant mass of diplioton can- 
didates after all selections for the combined 7 TeV and 8 TeV data 
sample. The inclusive sample is shown in (a) and a weighted version 
of the same sample in (c); the weights are explained in the text. The 
result of a fit to the data of the sum of a signal component fixed to 
nifi = 126.5 GeV and a background component described by a fourth- 
order Bernstein polynomial is superimposed. The residuals of the data 
and weighted data with respect to the respective fitted background 
component are displayed in (b) and (d). 



5.7. Results 

The distributions of the invariant mass, iriyy, of the 
diphoton events, summed over all categories, are shown 
in Fig. St a) and (b). The result of a fit including a signal 
component fixed to niH = 126.5 GeV and a background 
component described by a fourth-order Bernstein poly- 
nomial is superimposed. 

The statistical analysis of the data employs an un- 
binned likelihood function constructed from those of 
the ten categories of the 7 TeV and 8 TeV data samples. 
To demonstrate the sensitivity of this likelihood analy- 
sis. Fig. |5c) and (d) also show the mass spectrum ob- 
tained after weighting events with category-dependent 
factors reflecting the signal-to-background ratios. The 
weight Wi for events in category / e [1,10] for the 7 TeV 
and 8 TeV data samples is defined to be In (1 + 5,/B,), 
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where 5, is 90% of the expected signal for nin - 126.5 
GeV, and B, is the integral, in a window containing S 
of a background-only fit to the data. The values 5,7S, 
have only a mild dependence on nin- 

The statistical interpretation of the excess of events 
near niyy = 126.5 GeV in Fig. |4]is presented in Sec- 
tion |9] 

6. //-» WW**' -»ev//v channel 

The signature for this channel is two opposite-charge 
leptons with large transverse momentum and a large 
momentum imbalance in the event due to the escaping 
neutrinos. The dominant backgrounds are non-resonant 
WW, tt, and Wt production, all of which have real 
W pairs in the final state. Other important backgrounds 
include Drell-Yan events {pp^Zly^*^^ U) with Ef^"- 
that may arise from mismeasurement, W-i-jets events 
in which a jet produces an object reconstructed as the 
second electron or muon, and Wy events in which the 
photon undergoes a conversion. Boson pair produc- 
tion (Wy'/WZ**' and ZZ**') can also produce opposite- 
charge lepton pairs with additional leptons that are not 
detected. 

The analysis of the 8 TeV data presented here is fo- 
cused on the mass range 110 < mh < 200 GeV. It fol- 
lows the procedure used for the 7 TeV data, described 
in Ref. [106], except that more stringent criteria are ap- 
plied to reduce the W+jets background and some selec- 
tions have been modified to mitigate the impact of the 
higher instantaneous luminosity at the LHC in 2012. In 
particular, the higher luminosity results in a larger Drell- 
Yan background to the same-flavour final states, due to 
the deterioration of the missing transverse momentum 
resolution. For this reason, and the fact that the efi final 
state provides more than 85% of the sensitivity of the 
search, the same-flavour final states have not been used 
in the analysis described here. 

6.1. Event selection 

For the 8 TeV IVW**'— > evyuv search, the data 
are selected using inclusive single-muon and single- 
electron triggers. Both triggers require an isolated lep- 
ton with pf > 24 GeV. Quality criteria are applied 
to suppress non-collision backgrounds such as cosmic- 
ray muons, beam-related backgrounds, and noise in 
the calorimeters. The primary vertex selection fol- 
lows that described in Section |4] Candidates for the 
//— > WW**'— > ev/iv search are pre-selected by requir- 
ing exactly two opposite-charge leptons of diff'erent 
flavours, with px thresholds of 25 GeV for the leading 



lepton and 15 GeV for the sub-leading lepton. Events 
are classified into two exclusive lepton channels de- 
pending on the flavour of the leading lepton, where ep 
(pe) refers to events with a leading electron (muon). The 
dilepton invariant mass is required to be greater than 
10 GeV. 

The lepton selection and isolation have more strin- 
gent requirements than those used for the H — > ZZ**' — > 
A{ analysis (see Section |4]i, to reduce the larger back- 
ground from non-prompt leptons in the Mv final state. 
Electron candidates are selected using a combination of 
tracking and calorimetric information [85]; the criteria 
are optimised for background rejection, at the expense 
of some reduced efficiency. Muon candidates are re- 
stricted to those with matching MS and ID tracks il84l] , 
and therefore are reconstructed over I77I < 2.5. The 
isolation criteria require the scalar sums of the pi of 
charged particles and of calorimeter topological clus- 
ters within AR — 0.3 of the lepton direction (excluding 
the lepton itself) each to be less than 0. 12-0.20 times the 
lepton pj. The exact value differs between the criteria 
for tracks and calorimeter clusters, for both electrons 
and muons, and depends on the lepton pj. Jet selec- 
tions follow those described in Section |531 except that 
the JVF is required to be greater than 0.5. 

Since two neutrinos are present in the signal final 
state, events are required to have large £'™^\ E™'^^ is 
the negative vector sum of the transverse momenta of 
the reconstructed objects, including muons, electrons, 
photons, jets, and clusters of calorimeter cells not as- 
sociated with these objects. The quantity E™^^^ used 
in this analysis is required to be greater than 25 GeV 
and is defined as: ^i™!^^; = -E™"*^ sin A^min, where A^y^^ia 
is inin(A(^, |), and E"'^^ is the magnitude of the vec- 
tor E™^^ Here, A0 is the angle between E™^"^ and the 
transverse momentum of the nearest lepton or jet with 
Pi > 25 GeV. Compared to £'™'^\ rei increased 
rejection power for events in which the is™"* is gener- 
ated by a neutrino in a jet or the mismeasurement of an 
object, since in those events the E™^"* tends to point in 
the direction of the object. After the lepton isolation and 
^T rei requirements that define the pre-selected sample, 
the multijet background is negligible and the Drell-Yan 
background is much reduced. The Drell-Yan contribu- 
tion becomes very small after the topological selections, 
described below, are applied. 

The background rate and composition depend signif- 
icantly on the jet multiplicity, as does the signal topol- 
ogy. Without accompanying jets, the signal originates 
almost entirely from the ggF process and the back- 
ground is dominated by WW events. In contrast, when 
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produced in association with two or more jets, the signal 
contains a much larger contribution from the VBF pro- 
cess compared to the ggF process, and the background 
is dominated by it production. Therefore, to maximise 
the sensitivity to SM Higgs events, further selection cri- 
teria depending on the jet multiplicity are applied to the 
pre-selected sample. The data are subdivided into 0-jet, 
1-jet and 2-jet search channels according to the num- 
ber of jets in the final state, with the 2-jet channel also 
including higher jet multiplicities. 

Owing to spin correlations in the IVW**' system aris- 
ing from the spin-0 nature of the SM Higgs boson and 
the V-A structure of the W boson decay vertex, the 
charged leptons tend to emerge f rom the primary ver- 
tex pointing in the same direction llOTIl . This kinematic 
feature is exploited for all jet multiplicities by requiring 
that \lS.(pfc\ < 1-8, and the dilepton invariant mass, mu, 
be less than 50 GeV for the 0-jet and 1-jet channels. For 
the 2-jet channel, the mu upper bound is increased to 
80 GeV. 

In the 0-jet channel, the magnitude p^^ of the trans- 
verse momentum of the dilepton system, p^'^ - p^' +Pj^, 
is required to be greater than 30 GeV. This improves the 
rejection of the Drell-Yan background. 

In the 1-jet channel, backgrounds from top quark pro- 
duction are suppressed by rejecting events containing a 
i?-tagged jet, as determined using a Zj-tagging algorithm 
that uses a neural network and exploits the topology of 
weak decays of b- and c-hadrons lll08ll . The total trans- 
verse momentum, p^°^, defined as the magnitude of the 
vector sum p^' - p^' + p^^ + p^ + E™^\ is required 
to be smaller than 30 GeV to suppress top background 
events that have jets with pj below the threshold defined 
for jet counting. In order to reject the background from 
Z— > TT, the TT invariant mass, m^r, is computed under 
the assumptions that the reconstructed leptons are t lep- 
ton decay products. In addition the neutrinos produced 
in these decays are assumed to be the o nly s ource of 
£'™^'* and to be collinear with the leptons |ll09l. Events 
with \mrT - mz\ < 25 GeV are rejected if the collinear 
approximation yields a physical solution. 

The 2-jet selection follows the 1-jet selection de- 
scribed above, with the definition modified to in- 
clude all selected jets. Motivated by the VBF topol- 
ogy, several additional criteria are applied to the tag 
jets, defined as the two highest-px jets in the event. 
These are required to be separated in rapidity by a dis- 
tance \Ayjj\ > 3.8 and to have an invariant mass, mjj, 
larger than 500 GeV. Events with an additional jet with 
Pt > 20 GeV between the tag jets (yji < y < yji) are 
rejected. 



A transverse mass variable, nij is used to test 

for the presence of a signal for all jet multiplicities. This 
variable is defined as: 



)2- 



|p^^ -I- E™-^^p, 



where Ej 



p^p' 1^ -I- mjf. The statistical analysis of 
the data uses a fit to the mj distribution in the signal re- 
gion after the Acpu requirement (see Section l64] i. which 
results in increased sensitivity compared to the analysis 
described in Ref. lllllil . 

For a SM Higgs boson with nin = 125 GeV, the 
cross section times branching ratio to the evfiv final 
state is 88 fb for V* = 7TeV, increasing to 112 fb at 
= 8 TeV The combined acceptance times efficiency 
of the 8 TeV 0-jet and 1-jet selection relative to the ggF 
production cross section times branching ratio is about 
7.4%. The acceptance times efficiency of the 8 TeV 2-jet 
selection relative to the VBF production cross section 
times branching ratio is about 14%. Both of these fig- 
ures are based on the number of events selected before 
the final nij criterion is applied (as described in Sec- 
tion |631l. 

6.2. Background normalisation and control samples 

The leading backgrounds from SM processes produc- 
ing two isolated high-px leptons are WW and top (in 
this section, "top" background always includes both tt 
and single top, unless otherwise noted). These are es- 
timated using partially data-driven techniques based on 
normalising the MC predictions to the data in control 
regions dominated by the relevant background source. 
The W+jets background is estimated from data for all jet 
multiplicities. Only the small backgrounds from Drell- 
Yan and diboson processes other than WW, as well as 
the WW background for the 2-jet analysis, are estimated 
using MC simulation. 

The control and validation regions are defined by se- 
lections similar to those used for the signal region but 
with some criteria reversed or modified to obtain signal- 
depleted samples enriched in a particular background. 
The term "validation region" distinguishes these regions 
from the control regions that are used to directly nor- 
malise the backgrounds. Some control regions have sig- 
nificant contributions from backgrounds other than the 
targeted one, which introduces dependencies among the 
background estimates. These correlations are fully in- 
corporated in the fit to the nij distribution. In the fol- 
lowing sections, each background estimate is described 
after any others on which it depends. Hence, the largest 
background (WW) is described last. 
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6.2.7. W+jets background estimation 

The W+jets background contribution is estimated us- 
ing a control sample of events where one of the two lep- 
tons satisfies the identification and isolation criteria de- 
scribed in Section 16.11 and the other lepton fails these 
criteria but satisfies a loosened selection (denoted "anti- 
identified"). Otherwise, events in this sample are re- 
quired to pass all the signal selections. The dominant 
contribution to this sample comes from VK+jets events 
in which a jet produces an object that is reconstructed 
as a lepton. This object may be either a true electron or 
muon from the decay of a heavy quark, or else a product 
of the fragmentation identified as a lepton candidate. 

The contamination in the signal region is obtained by 
scaling the number of events in the data control sam- 
ple by a transfer factor. The transfer factor is defined 
here as the ratio of the number of identified lepton can- 
didates passing all selections to the number of anti- 
identified leptons. It is calculated as a function of the 
anti-identified lepton pj using a data sample dominated 
by QCD jet production (dijet sample) after subtracting 
the residual contributions from leptons produced by lep- 
tonic W and Z decays, as estimated from data. The 
small remaining lepton contamination, which includes 
Wy^*yWZ^*^ events, is subtracted using MC simulation. 

The processes producing the majority of same-charge 
dilepton events, W+jets, W/*VWZ**' and Z^*'Z^*', are 
all backgrounds in the opposite-charge signal region. 
W+jets and Wy'*' backgrounds are particularly impor- 
tant in a search optimised for a low Higgs boson mass 
hypothesis. Therefore, the normalisation and kinematic 
features of same-charge dilepton events are used to val- 
idate the predictions of these backgrounds. The pre- 
dicted number of same-charge events after the ^jj.g'j and 
zero-jet requirements is 216 + 7 (stat) + 42 (syst), while 
182 events are observed in the data. Satisfactory agree- 
ment between data and simulation is observed in vari- 
ous kinematic distributions, including those of A(pf( (see 
Fig. 13 a)) and the transverse mass. 

6.2.2. Top control sample 

In the 0-jet channel, the top quark background predic- 
tion is first normalised using events satisfying the pre- 
selection criteria described in Section ISTI This sample 
is selected without jet multiplicity or /j-tagging require- 
ments, and the majority of events contain top quarks. 
Non-top contributions are subtracted using predictions 
from simulation, except for W+jets, which is estimated 
using data. After this normalisation is performed, the 
fraction of events with zero jets that pass all selections 
is evaluated. This fraction is small (about 3%), since the 
top quark decay Wb has a branching ratio of nearly 
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Figure 5: Validation and control distributions for the 
//— » WW'*'— ► evjuv analysis. a) A0ff distribution in the same- 
charge validation region after the and zero-jet requirements, b) 
mx distribution in the WW control region for the 0-jet channel. The 
e/i and /le final states are combined. The hashed area indicates the 
total uncertainty on the background prediction. The expected signal 
for mu = 125 GeV is negligible and therefore not visible. 



1 . Predictions of this fraction from MC simulation are 
sensitive to theoretical uncertainties such as the mod- 
elling of initial- and final-state radiation, as well as ex- 
perimental uncertainties, especially that on the jet en- 
ergy scale. To reduce the impact of these uncertainties, 
the top quark background determination uses data from 
a ii-tagged control region in which the one-to-two jet ra- 
tio is compared to the MC simulation 1 1 12]. The result- 
ing correction factor to a purely MC -based background 
estimate after all selections amounts to 1 . 11 +0.06 (stat). 

In the 1-jet and 2-jet analyses, the top quark back- 
ground predictions are normalised to the data using con- 
trol samples defined by reversing the b-jet veto and re- 
moving the requirements on A(f)cc and mc(. The |Ayjj| 
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and mjj requirements are included in the definition of 
the 2-jet control region. The resulting samples are domi- 
nated by top quark events. The small contributions from 
other sources are taken into account using MC simula- 
tion and the data-driven W+jets estimate. Good agree- 
ment between data and MC simulation is observed for 
the total numbers of events and the shapes of the nij 
distributions. The resulting normalisation factors are 
1.11 ± 0.05 for the 1-jet control region and 1.01 + 0.26 
for the 2-jet control region. Only the statistical uncer- 
tainties are quoted. 

6.2.3. WW control sample 

The MC predictions of the WW background in the 
0-jet and 1-jet analyses, summed over lepton flavours, 
are normalised using control regions defined with the 
same selections as for the signal region except that the 
A(f>(( requirement is removed and the upper bound on 
m/f is replaced with a lower bound: m^f > 80 GeV. The 
numbers of events and the shape of the mx distribution 
in the control regions are in good agreement between 
data and MC, as shown in Fig. |5|b). WW production 
contributes about 70% of the events in the 0-jet control 
region and about 45% in the 1-jet region. Contamina- 
tions from sources other than WW are derived as for the 
signal region, including the data-driven W^-i-jets and top 
estimates. The resulting normalisation factors with their 
associated statistical uncertainties are 1.06 + 0.06 for the 
0-jet control region and 0.99 + 0.15 for the 1-jet control 
region. 

6.3. Systematic uncertainties 

The systematic uncertainties that have the largest im- 
pact on the sensitivity of the search are the theoretical 
uncertainties associated with the signal. These are de- 
scribed in Section |9] The main experimental uncertain- 
ties are associated with the JES, the jet energy resolu- 
tion (JER), pile-up, £'™^% the ^-tagging efficiency, the 
W+jets transfer factor, and the integrated luminosity. 
The largest uncertainties on the backgrounds include 
WW normalisation and modelling, top normalisation, 
and Wy^*'' normalisation. The 2-jet systematic uncer- 
tainties are dominated by the statistical uncertainties in 
the data and the MC simulation, and are therefore not 
discussed further. 

Variations of the jet energy scale within the system- 
atic uncertainties can cause events to migrate between 
the jet bins. The uncertainty on the JES varies from 
+2% to ±9% as a function of je t pj and rj for jets with 
Pt > 25 GeV and \t]\ < 4.5 llOlll . The largest impact of 
this uncertainty on the total signal (background) yield 



amounts to 7% (4%) in the 0-jet (1-jet) bin. The un- 
certainty on the JER is estimated from in situ measure- 
ments and it impacts mostly the 1-jet channel, where 
its effect on the total signal and background yields is 
4% and 2%, respectively. An additional contribution to 
the JES uncertainty arises from pile-up, and is estimated 
to vary between +1% and ±5% for multiple pp colli- 
sions in the same bunch crossing and up to ±10% for 
neighbouring bunch crossings. This uncertainty affects 
mainly the 1-jet channel, where its impact on the sig- 
nal and background yields is 4% and 2%, respectively. 
JES and lepton momentum scale uncertainties are prop- 
agated to the Ej"'^ measurement. Additional contri- 
butions to the Ej"^^ uncertainties arise from jets with 
Pj < 20 GeV and from low-energy calorimeter deposits 
not associated with reconstructed physics objects j 1 1 31] . 
The impact of the £'™^'* uncertainty on the total signal 
and background yields is ~3%. The efficiency of the b- 
tagging algorithm is calibrated using samples contain- 
ing muons reconstructed in the vicinity of jets II fill. 
The uncertainty on the b-jet tagging efficiency varies be- 
tween ±5% and ±18% as a function of the jet pj, and 
its impact on the total background yield is 10% for the 
1-jet channel. The uncertainty in the IV+jets transfer 
factor is dominated by differences in jet properties be- 
tween dijet and W+jets events as observed in MC sim- 
ulations. The total uncertainty on this background is 
approximately ±40%, resulting in an uncertainty on the 
total background yield of 5%. The uncertainty on the 
integrated luminosity is ±3.6%. 

A fit to the distribution of mj is performed in or- 
der to obtain the signal yield for each mass hypoth- 
esis (see Section 16.41 ). Most theoretical and experi- 
mental uncertainties do not produce statistically signif- 
icant changes to the mj distribution. The uncertainties 
that do produce significant changes of the distribution 
of mj have no appreciable effect on the final results, 
with the exception of those associated with the WW 
background. In this case, an uncertainty is included to 
take into account differences in the distribution of m j 
and normalisation observed between the MCFM 11 1511 . 
MC@NLO+HERWIG and POWHEG+PYTHIA gen- 
erators. The potential impact of interference between 
resonant ( Higg s-mediated) and non-resonant gg^ WW 
diagrams jl 1611 for mj > mn was investigated and found 
to be negligible. The effect of the WW normalisation, 
modelling, and shape systematics on the total back- 
ground yield is 9% for the 0-jet channel and 19% for 
the 1-jet channel. The uncertainty on the shape of the 
total background is dominated by the uncertainties on 
the normalisations of the individual backgrounds. The 
main uncertainties on the top background in the 0-jet 
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analysis include those associated with interference ef- 
fects between tt and single top, initial state an final state 
radiation, i>-tagging, and JER. The impact on the total 
background yield in the 0-jet bin is 3%. For the 1-jet 
analysis, the impact of the top background on the to- 
tal yield is 14%. Theoretical uncertainties on the Wy 
background normalisation are evaluate d for each jet bin 
using the procedure described in Ref. fll ITIl . They are 
+ 11% for the 0-jet bin and +50% for the 1-jet bin. For 
Wj* with nice < 7 GeV, a k-factor of 1 .3 + 0.3 is applied 
to the MadGraph LO prediction based on the compari- 
son with the MCFM NLO calculation. The k-factor for 
Wy*/WZ^*'> with nice > 7 GeV is 1.5 ± 0.5. These un- 
certainties affect mostly the 1-jet channel, where their 
impact on the total background yield is approximately 
4%. 



Table 5: The expected numbers of signal (rtiu = 125 GeV) and back- 
ground events after all selections, including a cut on the transverse 
mass of 0.15 mu < mj < m/y for mu = 125 GeV. The observed 
numbers of events in data are also displayed. The e/i and /le chan- 
nels ai'e combined. The uncertainties shown are the combination of 
the statistical and all systematic uncertainties, taking into account the 
constraints from control samples. For the 2-jet analysis, backgrounds 
with fewer than 0.01 expected events are marked with '-' . 
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0-jet 


1-jet 


2-jet 


Signal 


20 ±4 


5±2 


0.34 ±0.07 


WW 


101 ±13 


12±5 


0.10±0.14 


WZ^*"! /ZZ/Wy^*> 


12±3 


1.9±1.1 


0.10±0.10 


tt 


8±2 


6±2 


0.15±0.10 


tW/tb/tqb 


3.4±1.5 


3.7 ±1.6 




Z/y' +jets 


1.9± 1.3 


0.10±0.10 




W + jets 


15±7 


2± 1 




Total Background 


142 ±16 


26±6 


0.35 ±0.18 


Observed 


185 


38 






6.4. Results 

Table |5] shows the numbers of events expected from 
a SM Higgs boson with niH = 125 GeV and from the 
backgrounds, as well as the numbers of candidates ob- 
served in data, after application of all selection criteria 
plus an additional cut on mj of 0.75 m// < nij < nifj- 
The uncertainties shown in Table |5] include the system- 
atic uncertainties discussed in Section 16.31 constrained 
by the use of the control regions discussed in Sec- 
tion l6.2l An excess of events relative to the background 
expectation is observed in the data. 

Figure|6]shows the distribution of the transverse mass 
after all selection criteria in the 0-jet and 1-jet channels 
combined, and for both lepton channels together. 

The statistical analysis of the data employs a binned 
likelihood function constructed as the product of Pois- 
son probability terms for the efi channel and the jje 
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Figure 6: Distribution of the transverse mass, m-r, in the 0-jet and 
1-jet analyses with both e/i and /le channels combined, for events sat- 
isfying all selection criteria. The expected signal fornix = 125 GeV is 
shown stacked on top of the background prediction. The W+jets back- 
ground is estimated from data, and WW and top background MC pre- 
dictions are normalised to the data using control regions. The hashed 
area indicates the total uncertainty on the background prediction. 



channel. The mass-dependent cuts on nij described 
above are not used. Instead, the 0-jet (1-jet) signal re- 
gions are subdivided into five (three) mj bins. For the 
2-jet signal region, only the results integrated over mj 
are used, due to the small number of events in the final 
sample. The statistical interpretation of the observed 
excess of events is presented in Section |9] 



7. Statistical procedure 

The statistical proce dure used to interpret the data is 
described in Refs. llTlll 18l4l2lll . The parameter of in- 
terest is the global signal strength factor jj, which acts 
as a scale factor on the total number of events pre- 
dicted by the Standard Model for the Higgs boson sig- 
nal. This factor is defined such that fi - Q corresponds 
to the background-only hypothesis and /U = 1 corre- 
sponds to the SM Higgs boson signal in addition to the 
background. Hypothesised values of fi are tested with a 
statistic A(fj.) based on the profile likelihood ratio lll22ll . 
This test statistic extracts the information on the signal 
strength from a full likelihood fit to the data. The likeli- 
hood function includes all the parameters that describe 
the systematic uncertainties and their correlations. 

Ex clusi on limits are based on the CL, prescrip- 
tion lll23ll : a value of fi is regarded as excluded at 
95% CL when CL, is less than 5%. A SM Higgs bo- 
son with mass mn is considered excluded at 95% confi- 
dence level (CL) when fi - 1 is excluded at that mass. 
The significance of an excess in the data is first quan- 
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tified with the local po, the probability that the back- 
ground can produce a fluctuation greater than or equal 
to the excess observed in data. The equivalent formu- 
lation in terms of number of standard deviations, Z;, is 
referred to as the local significance. The global prob- 
ability for the most significant excess to be observed 
anywhere in a given searc h reg ion is estimated with the 
method described in Ref. u2M- The ratio of the global 
to the local probabilities, the trials factor used to correct 
for the "look elsewhere" effect, increases with the range 
of Higgs boson mass hypotheses considered, the mass 
resolutions of the channels involved in the combination, 
and the significance of the excess. 

The statistical tests are performed in steps of values 
of the hypothesised Higgs boson mass ;«//. The asymp- 
totic approximation |122] upon which the results are 
based has been validated with the method described in 
Ref. iH. 

The combination of individual search sub-channels 
for a specific Higgs boson decay, and the full combi- 
nation of all search channels, are based on the global 
signal strength factor fj. and on the identification of 
the nuisance parameters that correspond to the cor- 
related sources of systematic uncertainty described in 
Section!!] 



8. Correlated systematic uncertainties 

The individual search channels that enter the combi- 
nation are summarised in Table |6] 

The main uncorrected systematic uncertainties 
are described in Sections HHH for the //-> ZZ^*'^ 4f, 
yy and //— > WW**'— > ^v^v channels and in 
Ref. ilTIl for the other channels. They include the 
background normalisations or background model 
parameters from control regions or sidebands, the 
Monte Carlo simulation statistical uncertainties and 
the theoretical uncertainties affecting the background 
processes. 

The main sources of correlated systematic uncertain- 
ties are the following. 

1. Integrated luminosity: The uncertainty on the 
integrated luminosity is considered as fully correlated 
amon g channels and amounts to ±3.9% for the 7TeV 
data il32l [l33ll . except for the H^ZZ^*^^M and 
H^yy channels which were re-analysed; the uncer- 
tainty is + 1 .8% ll92ll for these channels. The uncertainty 
is +3.6% for the 8 TeV data. 

2. Electron and photon trigger identification: The 
uncertainties in the trigger and identification efficiencies 
are treated as fully correlated for electrons and photons. 



3. Electron and photon energy scales: The elec- 
tron and photon energy scales in the H^ZZ^*^^A( 
and yy channels are described by five parameters, 
which provide a detailed account of the sources of sys- 
tematic uncertainty. They are related to the calibration 
method, the presampler energy scale in the barrel and 
end-cap calorimeters, and the material description up- 
stream of the calorimeters. 

4. Muon reconstruction: The uncertainties affecting 
muons are separated into those related to the ID and MS, 
in order to obtain a better description of the correlated 
effects among channels using different muon identifica- 
tion criteria and different ranges of muon pj. 

5. Jet energy scale and missing transverse energy: 
The jet energy scale and jet energy resolution are af- 
fected by uncertainties which depend on the pi, rj, and 
flavour of the jet. A simplified scheme is used in which 
independent JES and JER nuisance parameters are asso- 
ciated with final states with significantly different kine- 
matic selections and sensitivity to scattering processes 
with different kinematic distributions or flavour com- 
position. This scheme includes a specific treatment 
for fe-jets. The sensitivity of the results to various as- 
sumptions about the correlation between these sources 
of uncertainty has been found to be negligible. An un- 
correlated component of the uncertainty on £'™'*'^ is in- 
cluded, in addition to the JES uncertainty, which is due 
to low energy jet activity not associated with recon- 
structed physics objects. 

6. Theory uncertainties: Correlated theoretical un- 
certainties affect mostly the signal predictions. The 
QCD scale uncertainties for ot//=125 GeV amount to 

for the ggF process, ± 1 % for the VBF and WH/ZH 
processes, and for the ttH process |[52i 153(1 : the 
small dependence of these uncertainties on m^ is taken 
into account. The uncertainties on the predicted branch- 
ing ratios amount to ±5%. The uncertainties related to 
the parton distribution functions amount to ±8% for the 
predominantly gluon-initiated ggF and ttH processes, 
and +4% for the predominantly quark-initiated VBF 
and WH/ZH processes l,7&. ,134-fl36 ]. The theoretical 
uncertainty associated with the exclusive Higgs boson 
production process with additional jets in the //— > yy, 
H-^ WW**'— > Mv and H — > r^r" channels is estimated 



using the prescription of Refs. 11531 ll 171 ll 1811 . with 



the noticeable difference that an explicit calculation of 
the gluon-fusion process at NLO using MCFM [I37I1 
in the 2-jet category reduces the uncertainty on this 
non-negligible contribution to 25 %. An additional 
theoretical uncertainty on the signal normalisation of 
+ 150%x(wH/TeV)^ (e.g. +4% for mn = 300 GeV) ac- 
counts for effects related to off-shell Higgs boson pro- 
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Table 6: Summary of the individual channels entering the combination. The transition points between separately optimised mu regions are indicated 
where applicable. In channels sensitive to associated production of the Higgs boson, V indicates a W or Z boson. The symbols igi and ffi represent 
direct products and sums over sets of selection requirements, respectively. 



Higgs Boson 
Decay 


Subsequent 
Decay 


Sub-Channels 


niH Range 
[GeV] 


J LA 

[fb-l] 


Ref. 


2011 V*=7TeV 


H ZZ^''> 


etvY 

e£qq 


{4e, 2e2/i, 2/j2e, 4p) 
[ee,^i^i] ^ {low, high pile-up| 
{fe-tagged, untagged) 


110-600 
200-280-600 
200-300-600 


4.8 
4.7 
4.7 


[87] 
[125] 
[126] 


H^yy 




10 categories {pxt <^Tfy^ conversion) ffi {2-jet) 


110-150 


4.8 


[127] 


H -» WW**' 


Mv 
tvqq' 


[ee, efiliie,)!)!] Ki (0-jet, 1-jet, 2-jet) 8i {low, high pile-up) 
{e,/j) (8 {0-jet, 1-jet, 2-jet) 


1 10-200-300-600 
300-600 


4.7 
4.7 


[106] 
[128] 


H ^TT 


Tlcp^had 
■'■had '"had 


{eti\ 8 {0-jet) ®[C(]<^{ 1-jet, 2-jet, VH] 
{e,/i) » {0-jet) igi {£™»" < 20 GeV, > 20 GeV) 
® {e,ii] igi (1-jet) ® [I] igi {2-jet) 
1 1-jet) 


110-150 
110-150 
110-150 


4.7 
4.7 
4.7 


[12^ 


VH -> Vhh 


Z — > VV 

W ^Cv 
Z->K 


£™'"' 6 { 120 - 160, 160 - 200, > 200 GeV) 
e {< 50, 50 - 100, 100 - 200, > 200 GeV) 
p\ e {< 50, 50 - 100, 100 - 200, > 200 GeV) 


110-130 
110-130 
110-130 


4.6 
4.7 
4.7 


[130] 


2012 V*^ =8 TeV 


H ZZ'-'^ 


At 


{4e, 2e2/i, 2^j2e, Ap.] 


110-600 


5.8 


[87] 


H^yy 




10 categories {pxt <^Tjy^ conversion) ffi {2-jet) 


110-150 


5.9 


[127] 


H WW**' 


evfjv 


\ep,pe\ 8 {0-jet, 1-jet, 2-jet) 


110-200 


5.8 


[131] 



Table 7: Characterisation of the excess in the H — » ZZ**' — » AC, //— > yy and //— > WW**'— > tvtv channels and the combination of all channels 
listed in Table|6] The mass value mi„ax for which the local significance is maximum, the maximum observed local significance Z/ and the expected 
local significance E(Zi) in the presence of a SM Higgs boson signal at mmn.v are given. The best fit value of the signal strength parameter fi at 
run = 126 GeV is shown with the total uncertainty. The expected and observed mass ranges excluded at 95% CL (99% CL, indicated by a *) are 
also given, for the combined = 7 TeV and = 8 TeV data. 



Search channel 


Dataset 


nimax [GeV] 


Z,[cr] 


E(Z,) [0-] 


fiiniH = 126 GeV) 


Expected exclusion [GeV] 


Observed exclusion [GeV] 




7 TeV 


125.0 


2.5 


1.6 


1.4± 1.1 






H ZZ<*' ^ A( 


8 TeV 


125.5 


2.6 


2.1 


1.1 ±0.8 








7 & 8TeV 


125.0 


3.6 


2.7 


1.2 ±0.6 


124-164, 176-500 


131-162, 170-460 




7 TeV 


126.0 


3.4 


1.6 


2.2 ±0.7 






H^yy 


8 TeV 


127.0 


3.2 


1.9 


1.5 ±0.6 








7 & 8TeV 


126.5 


4.5 


2.5 


1.8 ±0.5 


110-140 


112-123, 132-143 




7 TeV 


135.0 


1.1 


3.4 


0.5 ±0.6 






WW'*'^ Mv 


8 TeV 


120.0 


3.3 


1.0 


1.9 ±0.7 








7 & 8TeV 


125.0 


2.8 


2.3 


1.3 ±0.5 


124-233 


137-261 




7 TeV 


126.5 


3.6 


3.2 


1.2 ±0.4 






Combined 


8 TeV 


126.5 


4.9 


3.8 


1.5 ±0.4 






7 & 8TeV 


126.5 


6.0 


4.9 


1.4 ±0.3 


110-582 

1 13-532 (*) 


111-122, 131-559 
113-114, 117-121, 132-527 (*) 



duction and interference with other SM processes 115 311 . 

Sources of systematic uncertainty that affect both the 
7 TeV and the 8 TeV data are taken as fully correlated. 
The uncertainties on background estimates based on 
control samples in the data are considered uncorrected 
between the 7 TeV and 8 TeV data. 



9. Results 

The addition of the 8 TeV data for the ZZ^*'^ 4£, 
yy and //— > WW**'— > evyuv channels, as well as the 
improvements to the analyses of the 7 TeV data in the 



first two of these channels, bring a significant gain in 
sensitivity in the low-mass region with respect to the 
previous combined search 117,1 . 

9.1. Excluded mass regions 

The combined 95% CL exclusion limits on the pro- 
duction of the SM Higgs boson, expressed in terms of 
the signal strength parameter ji, are shown in Fig. |7ja) 
as a function of mn- The expected 95% CL exclu- 
sion region covers the nin range from 110 GeV to 
582 GeV. The observed 95% CL exclusion regions are 
11 1-122 GeV and 131-559GeV. Three mass regions 
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r(c) 



— Observed 

□ -2 In X{ii)<-\ 



m, [GeV] 

Figure 7: Combined searcli results: (a) The observed (solid) 95% CL 
limits on the signal strength as a function of m/y and the expec- 
tation (dashed) under the background-only hypothesis. The dark 
and light shaded bands show the ±lcr and ±2cr uncertainties on the 
background-only expectation, (b) The observed (solid) local po as a 
function of nifj and the expectation (dashed) for a SM Higgs boson 
signal hypothesis (/i = 1) at the given mass, (c) The best-fit signal 
strength ft as a function of ntfj- The band indicates the approximate 
68% CL interval around the fitted value. 



are excluded at 99% CL, 113-114, 117-121 and 132- 
527 GeV, while the expected exclusion range at 99% CL 
is 11 3-532 GeV. 

9.2. Observation of an excess of events 

An excess of events is observed near mH-l26 GeV in 
the ZZ'-*'>^ A{ and yy channels, both of which 
provide fully reconstructed candidates with high reso- 
lution in invariant mass, as shown in Figures Oa) and 
jS^b). These excesses are confirmed by the highly sen- 
sitive but low -resolution //— > WW^*'^ tviv channel, as 
shown in Fig.[8jc). 

The observed local po values from the combination 
of channels, using the asymptotic approximation, are 
shown as a function of nin in Fig.|7jb) for the full mass 
range and in Fig.|9]for the low mass range. 



- 2012 Exp. 

- 2012 Obs. 



2011-2012 Obs. 

= 7 TeV: jLdt = 
= 8 TeV: jLdt = 



4.7 fb"' 

5.8 fb' 



125 



130 



140 145 

niH [GeV] 

Figure 8: The observed local po as a function of the hypothesised 
Higgs boson mass for the (a) //— > ZZ'*'— ► 4f, (b) H~*yy and (c) 
H—>WW''^—>{v{v channels. The dashed curves show the expected 
local Po under the hypothesis of a SM Higgs boson signal at that mass. 
Results are shown separately for the = 7 TeV data (dark, blue), the 
V* = 8 TeV data (light, red), and their combination (black). 



The largest local significance for the combination of 
the 7 and 8 TeV data is found for a SM Higgs boson 
mass hypothesis of m//= 126.5 GeV, where it reaches 
6.0 cr, with an expected value in the presence of a SM 
Higgs boson signal at that mass of 4.9 cr (see also Ta- 
ble |7]i. For the 2012 data alone, the maximum lo- 
cal significance for the i/— » ZZ**' — > 4^, H—^yy and 
i/— > WW**'— » evfiv channels combined is 4.9 cr, and oc- 
curs at niH = 126.5 GeV (3.8 cr expected). 

The significance of the excess is mildly sensitive to 
uncertainties in the energy resolutions and energy scale 
systematic uncertainties for photons and electrons; the 
effect of the muon energy scale systematic uncertain- 
ties is negligible. The presenc e of t hese uncertainties, 
evaluated as described in Ref. lll38ll . reduces the local 
significance to 5.9 cr. 

The global significance of a local 5.9 cr excess any- 
where in the mass range 11 0-600 GeV is estimated to 
be approximately 5. 1 cr, increasing to 5.3 cr in the range 
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2011 -2012 _Qjj^ 
!_ is = 7 TeV: Jldt = 4.6-4.8 fb ' .... Exp. 

I- is = 8 TeV: Jldt = 5.8-5.9 fb ' □±1o 



1 

10"' 
10-2 
10"' 
lO" 
10 = 
10° 
10' 
10° 
10 = 
10'° 
10-" 

1 




10 115 120 125 130 135 140 145 150 

[GeV] 



Figure 9: The observed (solid) local po as a function of mu in the 
low mass range. The dashed curve shows the expected local po under 
the hypothesis of a SM Higgs boson signal at that mass with its ±1(T 
band. The horizontal dashed lines indicate the p- values corresponding 
to significances of 1 to 6 <t. 



110-150GeV, which is approximately the mass range 
not excluded at the 9 9% C L by the LHC combined SM 
Higgs boson search il39ll and the indirect constraints 
from t he g lobal fit to precision electroweak measure- 
ments ifiai . 

9.3. Characterising the excess 

The mass of the observed new particle is esti- 
mated using the profile likelihood ratio Aimn) for 
ZZ^*'>^ At and //-> yy, the two channels with the 
highest mass resolution. The signal strength is al- 
lowed to vary independently in the two channels, al- 
though the result is essentially unchanged when re- 
stricted to the SM hypothesis ji - 1. The leading 
sources of systematic uncertainty come from the elec- 
tron and photon energy scales and resolutions. The re- 
sulting estimate for the mass of the observed particle is 
126.0 + 0.4 (stat) + 0.4 (sys) GeV. 

The best-fit signal strength fi is shown in Fig. ^c) as 
a function of mn- The observed excess corresponds to 
fi - 1.4 ± 0.3 for niH - 126 GeV, which is consistent 
with the SM Higgs boson hypothesis ^ - \. A sum- 
mary of the individual and combined best-fit values of 
the strength parameter for a SM Higgs boson mass hy- 
pothesis of 126 GeV is shown in Fig. [TO] while more 
information about the three main channels is provided 
in Table H] 

In order to test which values of the strength and 
mass of a signal hypothesis are simultaneously consis- 
tent with the data, the profile likelihood ratio Aip., mn) is 
used. In the presence of a strong signal, it will produce 
closed contours around the best-fit point (/i, mn), while 



ATLAS 2011 - 2012 

W,Z H ^ bb 

Vs = 7TeV:|Ldl = 4.7 lb' 

H — > TT 

Vs = 7 TeV: jldl = 4.6-4.7 lb ' 

H ^ WW'"' ^ Iviv 

VE = 7TBV:|Ldl = 4.7lb' 
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H — > Yy 
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H = 1.4 ±0.3 



-1 1 

Signal strength ((x) 

Figure 10: Measurements of the signal strength parameter /i 
mfi=l26 GeV for the individual channels and their combination. 



for 



in the absence of a signal the contours will be upper 
limits on ju for all values of ot//. 

Asymptotically, the test statistic -2 In A(fj., nin) is dis- 
tributed as a distribution with two degrees of free- 
dom. The resulting 68% and 95% CL contours for the 
//— > yy and //— > WW*^*'— > Mv channels are shown in 
Fig -[HI where the asymptotic approximations have been 
validated with ensembles of pseudo-experiments. Sim- 
ilar contours for the //— » ZZ**'— > 4^ channel are also 
shown in Fig. [TT] although they are only approximate 
confidence intervals due to the smaller number of can- 
didates in this channel. These contours in the (fj., mn) 
plane take into account uncertainties in the energy scale 
and resolution. 

The probability for a single Higgs boson-like particle 
to produce resonant mass peaks in the //— >ZZ**'^4^ 
and //— > yy channels separated by more than the ob- 
served mass difference, allowing the signal strengths to 
vary independently, is about 8%. 

The contributions from the different production 
modes in the //— > yy channel have been studied in order 
to assess any tension between the data and the ratios of 
the production cross sections predicted in the Standard 
Model. A new signal strength parameter ju, is introduced 
for each production mode, defined by = cr,/cr, In 
order to determine the values of (/v,, ^j) that are simul- 
taneously consistent with the data, the profile likelihood 
ratio A{nj,Hj) is used with the measured mass treated as 
a nuisance parameter. 

Since there are four Higgs boson production modes at 
the LHC, two-dimensional contours require either some 
ju, to be fixed, or multiple ju, to be related in some way. 
Here, ji^gp and ii,jh have been grouped together as they 
scale with the ttH coupling in the SM, and are denoted 
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Figure II: Confidence intervals in the (p,mfj) plane for the 
H^ZZ'-'>^Ae, H^yy, and WW'^'>^Mv channels, including 
all systematic uncertainties. The markers indicate the maximum like- 
lihood estimates (fi, nifj) in the corresponding channels (the maximum 
likelihood estimates for ZZ^'>^ 4e and WW<*>^ Mv coin- 
cide). 



Figure 12: Likelihood contours for the H^yy channel in the 
0"ggF+rffl>j''VBF+v«) plane including the branching ratio factor 
B/BsM- The quantity /JjgF+n// (m^bf+vh) is a common scale factor 
for the ggF and tlH (VBF and VH) production cross sections. The 
best fit to the data {+) and 68% (full) and 95% (dashed) CL contours 
are also indicated, as well as the SM expectation (x). 



by the common parameter /iggp+z/T/- Similarly, jt/vBF and 
fivH have been grouped together as they scale with the 
WWH/ZZH coupling in the SM, and are denoted by the 
common parameter /ivBF+VH- Since the distribution of 
signal events among the 10 categories of the H^yy 
search is sensitive to these factors, constraints in the 
plane of figgP+nH x B/Bsm and //vbf+vh x B/Bsm, where 
B is the branching ratio for //— > yy, can be obtained 
(Fig.fTSI). Theoretical uncertainties are included so that 
the consistency with the SM expectation can be quanti- 
fied. The data are compatible with the SM expectation 
at the 1.5 cr level. 

10. Conclusion 

Searches for the Standard Model Higgs boson have 
been performed in the //— > ZZ***— > 4^, //— >yy and 
//-> WW***— > ev//v channels with the ATLAS experi- 
ment at the LHC using 5.8-5.9 fb ' of pp collision data 
recorded during April to June 2012 at a centre-of-mass 
energy of 8 TeV. These results are combined with ear- 
lier results fvf], which are based on an integrated lu- 
minosity of 4.6^.8 fb"' recorded in 2011 at a centre- 
of-mass energy of 7 TeV, except for the //— > ZZ**'— > 4£ 
and i/— > yy channels, which have been updated with the 
improved analyses presented here. 

The Standard Model Higgs boson is excluded at 
95% CL in the mass range 11 1-559 GeV, except for 
the narrow region 122-131 GeV. In this region, an ex- 
cess of events with significance 5.9 cr, corresponding 
to po - 1.7 X 10"^, is observed. The excess is driven 



by the two channels with the highest mass resolution, 
//— > ZZ***— > 4/" and H^yy, and the equally sensitive 
but low-resolution WW**'^ Mv channel. Taking 
into account the entire mass range of the search, 110- 
600 GeV, the global significance of the excess is 5.1 o", 
which corresponds to po = 1.7 x 10"^. 

These results provide conclusive evidence 
for the discovery of a new particle with mass 
126.0 ±0.4(stat) + 0.4 (sys) GeV. The signal 

strength parameter ju has the value 1.4 ± 0.3 at the 
fitted mass, which is consistent with the SM Higgs 
boson hypothesis fj. - I. The decays to pairs of vector 
bosons whose net electric charge is zero identify the 
new particle as a neutral boson. The observation in 
the diphoton channel disfavours the spin-1 hypothe- 
sis 1 14d 141 1. Although these results are compatible 
with the hypothesis that the new particle is the Standard 
Model Higgs boson, more data are needed to assess its 
nature in detail. 
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